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Abstract 

We calculate, using unitarity, a lower bound on the branching ratio 
XX ~^ 77 an d XX ~* 7-^> where x 1S an y halo dark matter particle that has 
W + W~ as one of the major annihilation modes. Examples of such par- 
ticles are supersymmetric particles with a dominant Higgsino component, 
or heavy triplet neutrinos. A substantial branching ratio is found for the 
77 and modes. We estimate the strength of the monoenergetic 7 ray 
lines that result from such annihilations in the Galactic or LMC halos. 



The question of the extent and nature of the dark matter in the Universe 
continues to be one of the most pressing problems of contemporary physics and 
astrophysics. Large scale structure and peculiar velocity observations seem to 
favour a value of the energy density much larger than what baryons can con- 
tribute due to nucleosynthesis constraints (for a recent review, see [0]). The 
recent candidate detections of microlensing events in the halo are, although 
very interesting, not conclusive as concerns the baryonic content of the dark halo 
of our galaxy (see, e.g., 0). Even if a large population of dark compact ob- 
jects (such as brown dwarfs) were formed at an early stage of galactic evolution, 
the overdensity of baryonic matter thus formed would accrete surrounding cold 
dark matter creating a particle dark matter halo density at least comparable in 
magnitude to the baryonic one. 

One of the favoured particle dark matter candidates is the lightest super- 
symmetric particle x> assumed to be a neutralino, i.e. a mixture of the super- 
symmetric partners of the photon, the Z° and the two neutral CP-even Higgs 
bosons present in the minimal extension of the supersymmetric standard model 
(see, e.g. 0]). The attractiveness of this candidate stems from the fact that 
its generic couplings and mass range naturally gives a relic density close to the 
critical one. Besides, its motivation from particle physics has recently become 
stronger due to the apparent need for 100 GeV - 10 TeV scale supersymmetry to 
achieve unification of the gauge couplings in view of recent LEP results B. 

When it comes to detecting cold dark matter particles, it seems that the tech- 
nology is not yet advanced enough for the ultrasensitive detectors developed to 
register the nuclear recoil and/or ionization to put interesting bounds on super- 
symmetric dark matter. On the other hand, indirect detection methods look quite 
promising. With large neutrino telescopes like DUMAND || and AMANDA [[7] 
now being deployed, there is a fair chance to detect energetic neutrinos from the 
center of the Sun or Earth coming from annihilations of captured neutralinos, if 
they constitute the dark matter halo ||. In fact, some bounds have already been 
obtained from the much smaller Kamiokande detector ||. 

The other way to detect particle dark matter in the halo is through indirect 
detection of positrons, antiprotons and 7 rays generated through the continuous 



annihilation of dark matter particles in the halo [10]. Observations of 7 rays have 



the advantage of giving a sensitive map of the galactic halo (since the annihilation 
rate depends on the square of the local dark matter density). In particular, if 
there is an enhancement of the density at the center of our galaxy [[11], [12], [13| 
or the Large Magellanic Cloud ||14j| , the 7 ray flux could stand out well above 
background. 

A particularly interesting annihilation process in the halo is XX ~ * 77 or 
XX ~ * Z 0/ ~f. Since these are two-body final states and the annihilating massive 
particles move with non-relativistic speed in the halo (typically v/c ~ 10 -3 ) the 



produced 7s will be nearly monoenergetic, meaning a 7 ray line signature [15 



The calculations of XX ~^ 77 f° r m x < m w have been made in quite some detail 
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[|T6| , [TT|j . The result is that the signal is only marginally detectable with present- 
day space detectors JL7|. On the other hand, it has recently been pointed out 
that a very massive \ in the TeV region could give a detectable signal in Air 



Cherenkov Telescopes (ACT) on the ground [fL8 |. These detectors have a very 
large effective area (on the order of 20 000 m 2 ) and a good proton rejection and 
energy resolution (~ 10 %) can be obtained with modern techniques. 

In fTf] the 77 process was estimated for the case of a nearly pure Bino, based 



on earlier calculations. However, the rates turn out to be very small [IS. 19]. 
There remains to do the more difficult calculation for a Higgsino, both for the 
77 and the Z°7 final states. The latter, which to our knowledge is considered 
for the first time here, will in fact turn out to be the most important one for 
Higgsinos (and for any other dark matter candidate that has a strong coupling to 
the W + W~ final state in its annihilations). Although a full calculation remains 
to be done for, say, the minimal supersymmetric extension of the standard model, 
we will use unitarity to put a strict lower bound on the 7 line signal as well as 
an estimate of the full rate. 

We consider first the Higgsino annihilation x°X° ~^ W + W~, mediated by 
chargino exchange in the t and u channel (there are also Z and H exchange 
contributions in general, but these vanish in the v — > limit). The calculation 
is identical for any Majorana particle that annihilates to W pairs through the 
exchange of a charged fermion. We write the x°X ± W T coupling (see p0| , ||J] for 
conventions; the index here identifies the lighter of the two charginos) 

r» = z {r(foP L + e p R ), (1) 

where P L>R = (1 ±7 5 )/2. 

In the limit where the lightest neutralino is a pure heavy Higgsino, the 
XX ~ * WW amplitude is dominated by the exchange of the lightest chargino, 
nearly degenerate with the neutralino, and 

|eo| = \fo\ = 1 (2) 
Projecting out the S wave part of the initial amplitude by the projector |^TJ 

we find the effective xxW + W~ vertex for polarization indices /x and v of the W 
bosons with four- momenta p + and 

V xx w + w- = { (e _° + /° 2) )e[^p-p + ], (4) 

V2V2(?7Co + m t+ - m w) ' 
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where e^ upa is the completeley antisymmetric constant tensor, and we use the 
abbreviation e[fiup_p + ] = e^ upa p f ' + p < L. This gives the annihilation rate 



12SiTm x mw (1 — u — Kf 
where u = (m x o/m^) 2 , k = (m x + / myy) 2 . This result coincides with the v — > 



limit of the expresion given in |20 



We are now prepared for the one-loop calculation of \X ~^ 77- We assume 
m x + >m x o >> m\y. It has been shown in a similar type of calculation [22] that 
it is a very good approximation, even when m x + ~ m x°-> ^° take ^ ne effective 
vertex (H) as pointlike : 



^xxW ~ Ce\pvk(k - qi - q 2 )}, (6) 

with 

r.-i 

2V2m 2 



(see Fig. 1). For the contribution of the imaginary part to the branching ratio 
which we will extract to obtain a lower bound on the cross section, expression (|6|) 
is even exact because the denominator of the chargino propagator is constant. 
One can easily convince oneself that the direct XXWW^f vertex (Fig. 1(b)) that by 
gauge invariance has to be present due to the derivative coupling in the effective 
vertex (|) does not contribute to our process. In addition, the ordinary contact 
WW'j'j term (Fig. 1(c)) contributes zero, since the intitial low- velocity XX state 
only projects out the part of the amplitude antisymmetric in the W momenta, and 
the contact term is symmetric in its indices and thus vanishes upon contraction. 

There remains to calculate the triangle graphs with Ws and in the customary 
linear gauge also Goldstone bosons circulating in the loop (Fig. 1(a)). Con- 
siderable simplification is obtained by choosing a non-linear gauge condition as 
originally suggested by Fujikawa |23|. This has been used in various applications 
such as calculating Higgs decays into two photons |24|, 25], the charge radius of 



the neutrino or radiative neutralino decay |]27||. The idea is to replace the 



usual gauge fixing term (we set £ = 1 for simplicity) 

C g .f. = ~\ {d^f - X - (d,Z» + m z G ) 2 - \8,W + » + im w G+\ 2 , (8) 

where G ±,CI are the Goldstone bosons eaten by and Z°, by 

= ~\ (W? - l - faz" + m z G ) 2 - + igW*)W+» + im w G + \\ (9) 

where is the neutral component of the SU (2) l gauge triplet (W^ = Z® cos 9 W — 
Afj,sin9 w ). The main advantage of this gauge is that the new contribution from 
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(||) to the trilinear ^W^G^ part of the Lagrangian actually cancels a correspond- 
ing piece of the original trilinear sector of the Lagrangian giving a vanishing 
total •yW ± G T coupling. In addition, the new terms entering the 7H^ ± H^ T and 
Z°W ± W T vertices, although making them superficially look more complicated 
(and less symmetric in the momenta of the three bosons) in reality make these 
vertices much simpler when at least one of the bosons is on mass shell. The 
diagram in Fig. 1(a) (and the one obtained by crossing the photon lines which 
actually gives an identical contribution) were calculated in the 't Hooft-Feynman 
version (£ = 1) of the nonlinear gauge (in fact, we also checked our results using 
the more cumbersome linear gauge). 

We find the amplitude for XX ~~ > 7(91)7(92) with photon polarization four 
vectors e M (<7i) and e u (q 2 ) to be 

_ -8ie 2 C r 4 m 2 x e[kfiu(q l + q 2 )] + k^kuq^} - k v t\k^q x q 2 \ 



(2tt) 4 J ^-m 2 w )((k-q 1 ) 2 -m 2 w )((k-q 1 -q 2 



rn 



w) 



where the ordinary prescription — > — ie is understood in the denomina- 
tor. 

Using standard techniques, this integral can be transferred to a symmetric 
integral in D = 4 — 2e dimensions after combining the factors in the denomi- 
nator using Feynman parametrization. The D-dimensional integrals are of two 
types, one convergent and one logarithmically divergent as D — > 4 (signalled by 
an 1/e divergence as Z) — > 4). The (infinite) renormalization constant is real 
and does not affect the imaginary part, which can be calculated in terms of ele- 
mentary functions since the 2-dimensional Feynman parameter integral becomes 
1-dimensional after using the 5-function coming from the ie term in the denom- 
inator. (This corresponds to using the Cutkosky rules to extract the imaginary 
part of the amplitude.) Writing 

Ce 2 

A7 = — ; — e[//z/gig 2 ] (M re + iM im ) , (11) 

we find 

M im = (3 2 log (i±| V (12) 



where (3 = Jl — {m w /m x ) 2 . 

In the limit of large \ masses (above, say, 500 GeV) this becomes 

M im ~ log ( . (13) 



in 



w . 



We have checked this result using still another method of calculation that relies 
on the generalized optical theorem, writing (the integral is over the WW phase 
space) 

M im = -J M XX ^ WW M;^ WW , (14) 
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where this time we employ the chargino propagator rather than the pointlike 
approximation for the xxWW vertex. This gives the same result and thus in 
addition proves the assertion that for the calculation of the imaginary part, the 
approximation of a pointlike effective vertex is exact whenever m x ± > m°. 

Once the imaginary (absorptive) part has been found, the real (dispersive) 
part can be calculated using dispersion relation techniques. However, for the 
pointlike vertex a subtraction will be needed, corresponding to a constant term 
resulting from the chargino propagator cutoff in the diagrams that generated 
the effective vertex. For very large x masses the leading logarithmic terms will 
dominate and can be calculated: 



= log(-i-), (15) 



m 



w 

where s = 4m?,. 

In the following we will use the imaginary part to get a strict lower bound for 
the rate estimates. To get an idea of the full result we will use the leading log 
estimate ( fTED for the real part. 

For the application we are considering, annihilation of heavy non-relativistic 
dark matter particles in the galactic halo, the 7 ray line coming from XX 77 
cannot realistically be discriminated from the corresponding line from XX ~^ Z ^. 
For example, the 7 energy for m x = 1 TeV is 1 TeV for the 77 mode and 
(4m 2 - m|)/(4m x ) = 998 GeV for the Z°7 mode, i.e. within the 10~ 3 - 10 -2 
spread of velocities of the annihilating particles (and of course well within the 
5 — 10 % energy resolution that at most can be achieved with present-day Air 
Cherenkov Telescopes). We should therefore add the 7 ray luminosities from the 
two sources. In fact, as we now shall see, the line strength from Z 0, j is expected 
to be a few times larger that that from the 77 final state. 

In the gauge we have chosen and to leading logarithmic accuracy, the difference 
between the 77 and Z 0/ y calculations is simply the replacement e 2 — > e 2 cot 9 W in 
Eq . (|TT|) . In the rate, this amounts to an enhancement by the factor cot 2 6 W ~ 3.4. 
When calculating the intensity of the line signal, the fact that each 77 event gives 
two photons is compensated by the division of the symmetry factor (= 2) due to 
the presence of two identical particles in the final state. 

In Fig. 2 we show our lower bound for the average number of 7 line photons 
of energy ~ m x per W + W~ annihilation: 

= av(xx W+W-) = Q - + COt 9w ) bg { Hi ) ' 



w 



as well as the value for this rate resulting from the adoption of (|16D for the 
dispersive part of the amplitude. When converting this into an effective line flux 
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per XX annihilation one has to correct by a factor giving the branching ratio of 
XX ~^ W + W~ normalized to the total annihilation cross section. This typically 
means a reduction of a factor of 2 for a Higgsino (due mainly to the ZZ channel), 
but no reduction for a species (e.g. a triplet neutrino [p8fl ) that does not couple 
to Z bosons at tree level. The remarkable feature of our result ( |iT|) is its very 
large size compared to the case of, e.g. a pure Bino [[18], As can be seen 
in Fig. 2, the effective branching ratio to line photons can easily be as large as 
10~ 2 . There will of course also be a continuous, diffuse flux of lower energy 7s 
coming from the fragmentation of the W + W~ and ZZ final states (see fl8|]). This 
distribution however lacks conspicuous features and has its main contribution at 
very low energies where the background 7 flux is much more severe and, what 
is worse, its shape is largely unknown. It may therefore prove very difficult to 
detect dark matter candidates through this "soft" photon flux. If detectors (e.g. 
ACTs) customized to achieve high energy resolution are built, a line signal would 
be an obvious feature to search for in the data. If such a line is found, there is 
certainly no known astrophysical background that could account for it, and the 
energy would simply correspond to the mass of the dark matter particle. 

Although we have focused on the Higgsino, one of the presently favoured 
dark matter candidates, our results are in fact of much greater generality for the 
following reasons. Any massive (TeV scale) neutral cold dark matter fermion 
may be expected to have a substantial annihilation strength into W + W~ pairs 
(although exceptions can be found, such as pure Binos). The fermion can be of 
Majorana or Dirac nature. For Majoranas, our analysis can be directly applied. 
The difference with Dirac fermions is that the CP properties do not forbid the 
slow fermions in the halo to annihilate in the triplet S (vector) state. However, 
due to gauge invariance, a massive spin 1 state can not decay into two photons, 
so the only contribution comes from the singlet S state which is the one we have 
already treated. Thus, we expect the logarithmic enhancement to be a quite 
general feature. The actual branching ratio line photons may, however, be lower 
by a factor of a few due to the new triplet contributions to the total annihilation 
rate. 

In the case of a heavy (m^ > rriyy) weak triplet Majorana neutrino (which 
was considered in J|8| for m N < m w ), our results need no modification, since 
WW is the only important tree level annihilation mode. That model is an inter- 
esting example of a viable particle dark matter model where direct detection is 
essentially impossible (since the neutrino does not have tree level Z couplings), 
but where 7 ray lines are predicted at a substantial level. 

We now turn to some estimates of the rates for the gamma line processes. 
Since the annihilation rate depends on the square of the neutralino number den- 
sity, the rate and angular distribution of the 7 line flux will depend on the model 
for the Galactic halo. In particular, there are some arguments that there may be 
a core region of the galaxy where the dark matter density is substantially higher 



than average In the model by Berezinsky et al |T2| the density profile falls 
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as 1/r 1,8 all the way from 5 Mpc down to the innermost part of the galaxy, where 
a black hole is thought to be residing. This model gives an enormous enhance- 
ment (several orders of magnitude) of the 7 ray flux from \ annihilations in the 
central region of the galaxy. Assuming this halo density profile, already existing 
measurements of high energy cosmic gamma ray fluxes constrain very strongly 
the supersymmetric models p9| (for a recent criticism of this model for the halo, 
see m. ' 



On the more conventional side, it has recently been pointed out [13] that the 
Large Magellanic Cloud may have its own halo (albeit with only partially known 
parameters), which possibly could give an enhanced flux from that direction. We 



will consider these cases, noting as in JT^] that for a smooth halo the rates we find 



are too small to be detectable with present technology. On the other hand, one 
of the the main advantages for turning the attention to the 7 ray signatures is 
just that this is a promising way of mapping the density distribution of the halo 
in the case it is not completely homogeneous. For instance, there is a possibility 
that string or texture generated perturbations in the cold dark matter density 
could have gone non-linear early and survived tidal disruptions to provide a very 



clumpy halo [[13 1 



For masses up to a few TeV, we can neglect galactic absorption of the 7 rays. 
According to |yj the optical depth at the 10 kpc length scale is around 10~ 3 . 



(Remember that for the high mass range we are considering, even for the Z7 
channel, E 1 = m x to an accuracy of within a percent.) 

For a localized source of dark matter annihilations at the center of the galaxy, 
the gamma line flux is given by 




where r ~ 8.5 kpc is the distance of the solar system from the center of the 
galaxy, and p(r) is the density profile within the core radius tq. For a generic dark 
matter particle that gives a substantial contribution to Q of the universe, (av ) to t 
is of the order of 10~ 26 cm 3 s _1 . Since the annihilation rate of neutralinos into 
W + W~ does not vanish nor become small when < v >— > 0, (cry ) to t is expected 
to be of the same order of magnitude in the galactic halo. As an example, this 
happens for a Higgsino of the minimal supersymmetric extension of the Standard 
Model of a mass around 1 TeV. Since we see from Fig. 2 that the number of line 
photons per annihilation is roughly of the order of a percent, we can use the value 
(av)n ne ~ 10~ 28 cm 3 s _1 in our estimates. 

In fact, for the case of a pure Higgsino, we can calculate the absolute value 
of the lower bound of the line cross section. In that case / = 1, eo = ±1 
in Eq. (1), and we find (assuming near Higgsino-chargino degeneracy, which 
is generally an excellent approximation for a nearly pure, massive Higgsino in 
minimal supersymmetric models) 
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7T(T 4 / 4m 2 \ 

2av{ XX - 77 ) + av( XX - Z^) > — log 2 , (19) 

x w V w / 

which corresponds to 0.6 • 10~ 28 cm 3 s _1 for m x = 1 TeV. Assuming, lacking a 
full calculation, a contribution from the dispersive part of the amplitudes of the 
same order of magnitude, we see that in this case the annihilation rate into line 
photons is indeed around 10 -28 cm 3 s _1 



Using the Ipser-Sikivie model for the core of the halo [ll|, a flux of around 
200 photons per year is predicted in an ACT of area 20 000 m 2 during a typical 
observation "year" of 2 • 10 6 s, with an estimated (but poorly known) background 
about an order of magnitude lower (including misidentified protons and electrons; 
see Jl8| for further discussions). This rate for Higgsino-like particles is two orders 



of magnitude larger than that expected for pure Binos [18|, T{|. In the model by 



Berezinsky et al [12|, the flux can be still at least two orders of magnitude larger. 



However, if one assumes a smooth halo distribution p(r) = p©(r 2 +a 2 )/(r 2 + a 2 ) 
with a ~ r Q , the rate falls below one event per year in our generic ACT detector 
in the direction of the galactic center. In the LMC model of Gondolo |14| the rate 
could be enhanced by an order of magnitude making it just about reaching the 
detectability limit. It should be noted, though, that the high threshold (200-300 
GeV) of present day ACTs make these processes rate limited. Going to larger 
areas and lower thresholds could greatly increase the discovery potential of this 
type of detector. 

To conclude, we have shown that the 77 and Z7 annihilation rates are quite 
high for any dark matter candidate that couples with full strength to W ± . A 
typical case is provided by a pure Higgsino, for which we have shown that the 7 
line strength is very much larger than that of a pure Bino. It should be noticed, 
though, that our calculations indicate this large 7 line rate to be generic for any 
dark matter candidate that couples to W ± with ordinary electroweak strength. If 
the galactic halo contains regions with higher dark matter density than the local 
(solar neighbourhood) value, the detection of high energy monoenergetic photons 
could at the same time determine the mass of the dark matter particle and map 
its galactic density distribution. 

This research has been sponsored by the Swedish Natural Science Research 
Council and a European Community Twinning grant. We wish to thank P. 
Gondolo, F. Martin and H. Rubinstein for useful discussions. We are especially 
grateful to M. Urban, whose idea to detect 7 lines from WIMP annihilation using 
Atmospheric Cerenkov Telescopes triggered this work. 



S 



References 

[1] See, e.g. M.S. Turner, Physica Scripta T36 (1991) 167. 

[2] C. Alcock et al, Nature 359 (1993) 393; E. Aubourg et al, Nature 365 (1993) 
623. 

[3] E. Gates and M.S. Turner, FERMILAB-PUB-93-357-A (1993). 

[4] H.E. Haber and G.L. Kane, Phys. Rep. 117 (1985) 75. 

[5] U. Amaldi, W. de Boer and H. Furstenau, Phys. Lett. B260 (1991) 447. 

[6] See, e.g., P.K.F. Greider, Europhysics News 23 (1992) 167. 

[7] See, e.g., D.M. Lowder et al., Nature 353 (1991) 331. 

[8] J. Silk, K. Olive and M. Srednicki, Phys. Rev. Lett. 55 (1985) 257; G.B. 
Gelmini, P. Gondolo and E. Roulet, Nucl. Phys. B351 (1991) 623; F. Halzen, 
M. Kamionkowski and T. Stelzer, Phys. Rev. D45 (1992) 4439. 

[9] M. Mori et al., Phys. Rev. D48 (1993) 5505. 

[10] J. Silk and M. Srednicki, Phys. Rev. Lett. 53 (1984) 624; F. Stecker, S. 
Rudaz and T. Walsh, Phys. Rev. Lett. 55 (1985) 2622; J. Ellis et al, Phys. 
Lett. B214 (1989) 403; M. Turner and F. Wilczek, Phys. Rev. D42 (1990) 
1001. 

[11] J.R. Ipser and P. Sikivie, Phys. Rev. D35 (1987) 3695. 

[12] V.S. Berezinsky, A.V. Gurevich and K.P Zybin, Phys. Lett. B294 (1992) 
221. 

[13] J. Silk and A. Stebbins, Astrophysical Journal 411 (1993) 439. 

[14] P. Gondolo, contributed paper to TAUP 93, Gran Sasso, Italy; to appear in 
the Proceedings. 

[15] L. Bergstrom and H. Snellman, Phys. Rev. D37 (1988) 3737; S. Rudaz, Phys. 
Rev. D39 (1989) 3549; A. Bouquet, P. Salati and J. Silk, Phys. Rev. D40 
(1989) 3168. 

[16] R. Barbieri and G. F. Giudice, Nucl. Phys. B306 (1988) 63; L. Bergstrdm, 
Phys. Lett. 225B (1989) 372. 

[17] V.S. Berezinsky, A. Bottino and V. de Alfaro, Phys. Lett. B274 (1992) 122. 

[18] M. Urban, A. Bouquet, B. Degrange, P. Fleury, J. Kaplan, A.L. Melchior 
and E. Pare, Phys. Lett. B293 (1992) 149. 



9 



[19] Note that there is an error in Eq. (20) of |]I8 |, entirely due to one of the 
present authors (JK). The prefactor should read 8 x 10 -31 em's -1 instead of 
4.2 x 10~ 29 em's" 1 . 

[20] K. Griest, M. Kamionkowski and M. Turner, Phys. Rev. D41 (1990) 3565. 
[21] J.H. Kiihn, J. Kaplan and O. Safiani, Nucl. Phys. B157 (1979) 125. 
[22] L. Bergstrdm, Phys. Lett. 225B (1989) 372. 
[23] K. Fujikawa, Phys. Rev. D7 (1973) 393. 

[24] M.B. Gavela, G. Girardi, C. Malleville and P. Sorba, Nucl. Phys. B193 (1981) 
257. 

[25] L. Begstrdm and G. Hulth, Nucl. Phys. B259 (1985) 137. 

[26] N.M. Monyonko and J.H. Reid, Prog. Theor. Phys. 73 (1984) 734. 

[27] H.E. Haber and D. Wyler, Nucl. Phys. B323 (1989) 267. 

[28] P. Chardonnet and P. Salati, Phys. Lett. B262 (1991) 307. 

[29] M. Urban, A. Bouquet, B. Degrange, M. Chantell, P. Fleury, J. Kaplan, A.L. 
Melchior, E. Pare, X. Sarazin and T.Weekes, Proceedings of the International 
Conference on Cosmic Ray Physics (Calgary, July 1993), preprint PAR- 
LPTHE 93-19; V.S. Berezinsky, A. Bottino and G. Mignola, Gran Sasso 
preprint LNGS 94/90 (1994). 

[30] R.A. Flores and J. Primack, UC Santa Cruz preprint SCIPP 93/01 (1994). 

[31] V.S. Berezinsky et al., Astrophysics of Cosmic Rays, North-Holland, 1990. 



10 



gure Captions 

1. Some diagrams relevant to the annihilation XX ~^ 77- F° r slow Majorana 
particle annihilation, only the diagram in (a) contributes (the diagrams 
obtained by crosing the photon lines are not shown in the figure). 

2. The average number F 7 (Eq. fllTD ) of 7 line photons from XX ~~ * 77 an d Z°7 
normalized to the W + W~ annihilation rate, as a function of the x mass. 
The solid line is the lower bound obtained by using only the imaginary part 
of the amplitude which is given by unitarity. The dashed line is an estimate 
of the full rate obtained by taking the leading logarithmic result of Eq. (|16|) 
for the real (dispersive) part of the amplitude. 
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